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EXTRACTION FROM FLIGHT DATA OF LONGITUDINAL AERODYNAMIC 
COEFFICIENTS IN MANEUVERING FLIGHT FOR F-8C AIRCRAFT 


William T. Suit 
Langley Research Center 

SUMMARY 

Flight -test data have been used to extract the longitudinal aerodynamic parameters 
of the F-8C aircraft at moderate to high angles of attack. The data obtained were from 
perturbations of the aircraft from steady turns , with trim normal accelerations from 
1.5g to 3.0g. The angle- of-attack variation from trim was approximately ±4° and maxi- 
mum angle of attack reached during the tests was approximately 16°. 

Although wind-tunnel data indicate that the lift and pitching moments are somewhat 
nonlinear with angle of attack, the linear aerodynamics extracted from the flight tests did 
permit computation of motion time histories which were in close agreement with the 
measured time histories. The aerodynamic parameters extracted from flight data were 
in reasonably good agreement with other data available for a lg condition. 

INTRODUCTION 

The National Aeronautics and Space Administration is currently involved in research 
on fly-by-wire control systems for aircraft. A fly-by-wire system replaces a direct-link 
mechanical control system with an electronic system to activate controls. Signals from 
the pilot go to a digital control system which commands the positions of the controls. 

These commands may be modified by a stability augmentation system. However, to design 
an augmentation system, the aerodynamic parameters of the aircraft must be known. The 
aircraft presently used by NASA to study digital fly -by -wire systems is an F-8C aircraft. 
The F-8C used as a test bed for digital fly -by -wire systems had a standard airframe. The 
aerodynamics available for the subject aircraft were determined primarily from wind- 
tunnel and analytical results and the mathematical model of the aircraft was considered to 
be reasonable, especially for trimmed level flight. To substantiate the existing model for 
the F-8C at moderate to high angles of attack, some flight tests were made with the angle 
of attack as high as 16°. 

A maximum -likelihood extraction procedure was used to examine the flight data. In 
this procedure, a set of equations of motion is used to calculate aircraft response to 



specified control inputs. Initial estimates of aerodynamic parameters (from theory or 
from wind-tunnel tests) are used for the initial motion computations. The maximum- 
likelihood extraction program then iterates on ithe aerodynamic parameters to select a 
set that maximizes a conditional likelihood function. This program has been used to 
determine the aerodynamic parameters for several aircraft in the lg trimmed flight 
condition. (See refs. 1 to 3.) The details of the program are contained in reference 4. 
The program has not been used previously at the Langley Research Center for flight 
data taken at moderate to high angles of attack where the trimmed load factor is greater 
than lg. 

The purpose of this paper is to present the longitudinal aerodynamic parameters of 
the F-8C aircraft as obtained from flight data flown at Mach numbers 0.7 and 0.8 and in 
trim with a load factor of 1.5g to 3.0g. The equations used and additional information on 
the flight data are presented and are followed by some comments on the extraction proce- 
dure and a discussion of the results of the study. 

SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 

a acceleration, g units 

c wing mean geometric chord, m (ft) 

g acceleration due to gravity, g units 

h altitude, m (ft) 

I moment of inertia, kg-m^ (slug-ft^) 

l f. distance from aircraft center of gravity to mean aerodynamic chord of hori- 

zontal tail, m (ft) 

M Mach number 

P period of oscillatory motion, sec 

p rate of roll, rad/sec or deg/sec 


2 


i in i 


ii i min ii 


i i 


q 


rate of pitch, rad/ sec or deg/ sec 


r rate of yaw, rad/sec or deg/sec 

5 wing area, (ft^) 

tjy '2 time to damp to half -amplitude, sec 

u component of velocity along X body axis, m/sec (ft/sec) 

V aircraft total velocity, m/sec (ft/sec) 

v component of velocity along Y body axis, m/sec (ft/sec) 

Wj. aircraft weight, N (lb) 

w component of velocity aong Z body axis, m/sec (ft/sec) 

ith component of state vector X 
a angle of attack, rad or deg 

6 e tail-plane deflection, positive for trailing edge down, rad or deg 

£ damping ratio 

6 pitch angle, rad or deg 

p air density, kg/m^ (slugs/ft^) 

cp bank angle, rad 

Coefficients: 

C m pitching-moment coefficient 

axial-force coefficient, positive along X body axis direction 
Cg normal-force coefficient, positive along Z body axis direction 
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c„ = 

8C Z 

c - 8C m 


da 

Cm 6 e " 86 e 

c m = 

U1 a 

8C m 

30 ? 

dCy 
Cy - f 

z< i a at 
2V 

Cv = 

9C X 

Q _ 8C Z 


da 

Z 6 e 95 e 

r> 

dC m 

r> _ 9 ^m 

U mq - 

2V 

2V 


Subscripts: 
c computed 

e tail plane 

f measured flight 

o indicates coefficient at trim conditions 

t indicates state at trim conditions 

X,Y,Z body -coordinate axes 

A dot over a symbol signifies a derivative with respect to time. 
Symbols used on computer plots: 

ALPHA = a - rad 

AXI acceleration along X body axis, g units 

AZI acceleration along Z body axis, g units 

DE = 6 e " 5 e,t> rad 

G acceleration due to gravity, 32.2 ft/sec2 (9.8 m/sec2) 
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Q 

rate of pitch, 

rad/sec 


THETA 

pitch angle, rad 



U 

velocity along 

X 

body axis, m/sec 

(ft/sec) 

V 

velocity along 

Y 

body axis, m/sec 

(ft/sec) 

w 

velocity along 

Z 

body axis, m/sec 

(ft/sec) 




EQUATIONS OF MOTION 


The equations of motion used in this study are referred to a body-axis system shown 
in figure 1 and are as follows: 

X-force: 

u = rv - qw - g sin 9 + |pV 2 S ^-[c x>0 + C X a ( a ~ “t) 

Z -force: 


w = qu - pv + g cos 9 cos cp + ~ pV 2 S -§~ 

l Wt 


C Z,o + C Z Q! (“" “t) 


+ C Z 6 ( 6 e - 6 e,t) + c Z q lv 


Pitching moment: 


4 = 


lz ,' Ix pr + ^(r2 - p2) ♦ 
r Y L Y 


ip Y?sc 
2 P I v 


C 


m,o 



t) 


+ C m „ + 

m q 2 y 


'm 


ac 
a 2V 


+ C 



a x = |(u + qw - rv + g sin 9 ) 
a z = |(w + pv - qu - g cos 9 cos cp) 
9 = q cos cp - r sin cp 
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FLIGHT DATA 


Description of Airplane 

The subject aircraft is a modified prototype F-8C and has been a flight -test vehicle 
since its manufacture in 1958. This aircraft is a single-seat high-performance airplane 
with a single jet engine embedded in the fuselage. Pitch control is achieved with a unit 
horizontal tail. The center of gravity used was at 29.0-percent mean geometric chord. 

The X body axis was parallel to and 10.16 cm (4 in.) above water line 100. (See fig. 1.) 
The geometric characteristics of the aircraft are given in table I. 

Flight Tests 

The data used in this report were obtained from flights flown at the NASA Flight 
Research Center. The pilot was instructed to fly a coordinated turn at nominal Mach 
numbers of 0.7 and 0.8 with nominal trim angles of attack of 9° and 13°. As can be seen 
from figures 2 to 5, the sideslip angle was less than ±2° at all times and the aircraft 
returned to the initial trim angle of attack after being perturbed by elevator steps and 
pulses. The angle-of-attack variation from trim was approximately ±4°. The aircraft 
stability augmentation systems were off during the tests. The actual test conditions for 
each individual run are given in table II. The aircraft weights and inertias listed in 
table II were obtained from tables supplied by NASA Flight Research Center and these 
weights and inertias were calculated as a function of the percent fuel in the aircraft. The 
weights and inertias used are average values for the test duration. Since the weight and 
inertias varied less than 3 percent and 1 percent, respectively, these variations were 
not accounted for in the parameter estimation. 

The data pertinent to this study, which were recorded during the flight tests, included 
longitudinal acceleration, a.yj normal acceleration, a n ; the difference between total pres- 
sure and static pressure as measured on a nose boom extending 1.83 m (6 ft) in front 
of the airplane; pitch attitude, 9; bank angle, cp; pitch rate, q; yaw rate, r; roll rate, p; 
indicated angle of attack, a and indicated angle of sideslip, /3, as measured by vanes on 
the nose boom; pressure altitude; control surface positions (aileron 5 a , elevator 6 e , and 
rudder 6 r ); and time, t. The full-scale range of the instruments is given in table III. 

All the data were stored on an onboard magnetic tape recorder by using a pulse control 
modulation (PCM) recording system. Additional information on the data acquisition sys- 
tem can be found in reference 4. 


Data Preparation 

The data were initially recorded, digitized, and calibrated at NASA Flight Research 
Center. A digital tape with the data in engineering units was sent to Langley Research 
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Center. The accelerometers were corrected for instrument location. The difference in 
the measured total pressure and static pressure was assumed to be the dynamic pressure. 
Density was determined from the standard atmosphere tables for the measured pressure 
altitude and the airspeed was calculated. The indicated a and /3 were corrected for 
the effects of aircraft angular rates by using the 0.29c center of gravity. The linear 
velocities along the vehicle body axes were calculated from the airspeed, angle of attack, 
and angle of sideslip. All the data were put on a tape at a rate of 20 points per second. 
The tape was then ready for use in the extraction program. 

PARAMETER-ESTIMATION PROCEDURE 


The parameter-estimation procedure used in this study is an iterative procedure 
which maximizes the conditional likelihood function L (aerodynamic parameters, 
weights, and initial conditions): 

N 

I ( x i,f - x i,c) 
i= 1 


(2tt) 1/2 |R| 1/2 


exp 


where R is the estimate of the error covariance matrix and X is the vector describing 
the state of the aircraft. The states used in the likelihood function were u, w, q, 0, 
a z , and a x . The calculated states c were determined by using the equations of 
motion previously introduced. In these equations the lateral quantities p, r, v, and cp 
were input directly into the equations from the flight data tape. The weighting matrix which 
is R - 1 can be obtained from the complete error covariance matrix, the diagonal terms of 
the error covariance matrix, or a diagonal matrix with fixed weights on the diagonal, at the 
discretion of the investigator. The diagonal form was used in this investigation and there- 
fore, the weights represent the estimated lower bound of the noise on the measured slates. 
The use of the likelihood function in parameter identification is discussed in reference 5. 
Maximizing the likelihood function results in a parameter-updated equation which is given 
by 


AP = 


(m t r _1 m) 


-1 


M T X; 


where M is the matrix of sensitivities of the calculated states with respect to the 
unknown parameters, T denotes the transpose and -1 denotes the inverse. (See ref. 5.) 
The matrix is the estimated parameter covariance matrix. The updated equation is deter- 
mined by forming a set of differential equations with the changes in the unknown param- 
eters as the variables. This set of simultaneous equations is then solved by the least- 
squares method to give the updated equation. (See ref. 5.) 
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After the convergence of the likelihood function, for a given flight data record, the 
extracted aerodynamic derivatives were examined. The derivatives were accepted as 
well determined if (a) the standard deviations of the computed time histories of the air- 
craft motion from the measured time histories were less than 3 percent of the full-scale 
range of the instrument used to measure the quantity examined, (b) the changes in the 
derivatives were less than one-hundredth of the derivative value for successive iterations, 
and (c) the estimated standard deviation of each derivative was less than about one-tenth 
of the extracted value of the derivative. 


Past experience has indicated a very high correlation between C m and C m • . 

4 ^ 

To eliminate this correlation, a value for Cm^ was chosen that was similar to values 

of C m . that were determined in reference 6. The value chosen was -0.6 and C m . was 
1X1 a 11 a 

held fixed at this value for all the runs. Also the geometric relation C 7c = C„, £- 

6e 6e It 

was assumed and so C 7 was not extracted but was calculated from the extracted value 

Z 8e 


RESULTS AND DISCUSSION 

Data for the flight conditions listed in table II were used with the model given to 
determine iteratively a set of aerodynamic derivatives for each of the flight conditions. 

The measured and computed time histories for each flight condition are shown in figures 2 
to 5. The measured data are represented by dotted lines in the figures. The computed 
time histories shown are those attained after the differences between the measured and 
calculated trajectories become constant. The figures show that in all cases the com- 
puted time histories were generally in close agreement with the flight records. Table IV 
gives the standard deviations of the computed states from the measured states. The 
inverses of the quantities in table IV were the diagonal terms used in the weighting matrix 
when the fits to the flight data shown in figures 2 to 5 were obtained. The standard devia- 
tions of the individual fits can be seen to be less than 3 percent of full scale, which was 
the uncertainty in the measured data (see table III); and in many cases, the standard 
deviations were less than 1 percent of the full-scale measured quantity. 

The derivatives extracted for each flight condition (the derivatives which resulted in 
the computed time histories of figs. 2 to 5) are listed in table V along with their estimated 
standard deviations. The correlation matrices for the aerodynamic parameters extracted 
are shown in table VI. An examination of table V indicates that the numbers extracted 
seem to be reasonable. The estimated variances of the parameters were less than 10 per- 
cent of the values extracted; thus confidence in the values obtained was indicated. 

For comparison, values of some aerodynamic coefficients obtained from reference 6 
are also shown in table V. The numbers given are for an altitude of 12.19 km (40 000 ft) 
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and were transformed from stability to body axes and converted to a 0.29c center-of- 
gravity location. 

The wind-tunnel and analytical investigation used for comparison showed a some- 
what nonlinear variation with pitching moment, drag force, and lift force with angle of 
attack for the angle -of-attack range of this investigation. These nonlinearities were not 
considered severe enough to make a linear aerodynamic model inappropriate. To show 
this, the curves of pitching moment, drag, force, and lift force plotted against angle of 
attack were obtained from reference 6 and straight-line approximations were made to 
these curves at the trim angles of attack for the flight tests. The derivatives estimated 
by use of these approximate curves were in reasonable agreement with the extracted 
derivative values. 

The results of this investigation seem to substantiate the curves generated from 
analytical and wind-tunnel studies for the test angle- of -attack range. The results also 
show that for the flight data examined, a mathematical model using linear aerodynamics 
was adequate to describe the response motions. To determine how well the numbers from 
reference 6 would represent the aircraft in the flight conditions of the subject study, the 
transformed numbers were put into the mathematical model and used to fit the two runs 
with the normal accelerations closest to lg. The fits to the data are shown in figures 6 
and 7. As can be seen, the fits to the data are not bad, but are not as good as those for 
the mathematical model determined by the extraction process. 

The value of C mQ from reference 6 was different from the value extracted; how- 
ever, a corresponding difference in the damping between the calculated response using the 
parameters from reference 6 and the calculated response of the extracted model was not 
observed. The flight data and extracted model appeared to have about 10 percent more 
damping than the calculated response using a model based on the parameters of refer- 
ence 6. In an attempt to get a better comparison of the two models and to substantiate the 
extracted model, the period and time to damp to half -amplitude of the short-period mode 
were determined by examining the damping envelope of the flight data and also calculating 
the period P and time to damp to half -amplitude ^\/2 ky use of the values from table V 
and equations from reference 7. The values of t, P, and £ are given in table VII. As 
can be seen the stability characteristics of the aircraft, as defined by the period and time to 
damp to half-amplitude for the longitudinal short period mode, are similar for the extracted 
parameters and for the parameters from reference 6. In this case a simplified analysis has 
shown that even though the parameter values in the two models, and especially C^, are 
different the response motions they define are similar. 

Additional verification of the aerodynamic parameters of table V was obtained by 
using, them to calculate motions , and comparing these calculated motions with motions 
from additional runs flown at the same time as the test runs but not used during the 
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extraction process. These comparisons are shown in figures 8 and 9. It can be seen 
that the fits to the additional data are essentially as good as the fits to the data used when 
parameters were extracted. 


CONCLUDING REMARKS 

Flight-test data flown at high to moderate angles of attack and using the F-8C air- 
plane has been used to determine the longitudinal aerodynamic parameters of the aircraft 
at four flight conditions. The tests were conducted with the aircraft trimmed in a steady 
turn with angles of attack of approximately 9° and 13° and at Mach numbers of 0.7 and 0.8. 

The extracted parameters resulted in a good fit to the flight data and the parameter 
values obtained were in fair agreement with values obtained from wind-tunnel and analyt- 
ical studies. The extracted aerodynamic model was further verified by comparing the 
period and time to damp to half -amplitude as calculated by using the extracted parameter 
values with the period and time to damp to half-amplitude measured from the flight -data 
traces. The extracted mathematical model was also used to fit data not used during the 
extraction process and the resulting fit was essentially as good as the fit obtained when 
the parameters were extracted. 

The wind-tunnel and analytical investigations indicated that the nonlinear variation 
of pitching moment and lift force with angle of attack for the angle-of-attack range of the 
investigation was not severe enough to make a linear aerodynamic model inappropriate. 
The results seem to substantiate this conclusion and show that for the flight data exam- 
ined, a mathematical model using linear aerodynamics was adequate to describe the 
response motions at the test angles of attack. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 10, 1975 
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TABLE L- GEOMETRIC CHARACTERISTICS OF THE F-8C AIRCRAFT 


Fuselage: 

Length, m (ft) 16.52 (54.17) 

Wing: 

Area, m 2 (ft 2 ) 34.88 (375) 

Aspect ratio 3.4 

Span, m (ft) 10.88 (35.67) 

Mean geometric chord, m (ft) 3.59 (11.78) 

Vertical tail: 

Area, m 2 (ft 2 ) 10.14 (109) 

Aspect ratio 1.5 

Span, m (ft) 3.89 (12.75) 

Rudder : 

Area, m 2 (ft 2 ) 1.17 (12.56) 

Horizontal tail: 

Area, m 2 (ft 2 ) 8.68 (93.4) 

Aspect ratio 3.5 

Span, m (ft) 5.52 (18.1) 

Tail length, center of gravity to quarter- chord point 

of mean geometric chord, m (ft) 5.06 (16.6) 
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TABLE n.- FLIGHT CONDITIONS 


Nominal altitude 

Nominal 

Mach 

number 

Nominal 

“t* 

deg 

Trim 

elevator, 

Trim bank 
angle, 

m 

ft 

6 e,t> 

deg 

deg 

10 370 

34 000 

0.71 

9.2 

-8.13 

50 

10 370 

34 000 

.66 

13.0 

-9.74 

60 

10 370 

34 000 

.81 

8.2 

-7.45 

63 

10 370 

34 000 

.78 

12.0 

-10.3 

70 







Inertias 




weigni 

ix 

I 

Y 

iz 

: XZ 

N 

lb 

kg-m 2 

slug- ft 2 

kg-m 2 

slug- ft 2 

kg-m 2 

slug- ft 2 

O 

kg-m 

slug- ft 2 

*94 000 

21 125 

12 500 

9200 

118 000 

86 800 

124 000 

91 600 

4030 

2970 


*Any errors in assuming nominal conditions were no greater than 3-percent sys- 
tem uncertainty on the measurements themselves. 
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TABLE m.- INSTRUMENT RANGES 

u and w are calculated from u = V cos a and w = V sin a. 
Individual sensors are basically more accurate than 3 percent 
of full scale; however, because of unknown errors, the effects 
of incompatibilities between measured states and processing 
errors the system accuracy was assumed to be 3 percent of the 
full-scale range of the instrument for the data as used during 
the extraction procedure. 


State 

Range 

a X 

±0.5g 

a Z 

-4g to +6g 

V 

30.91 to 515.15 m/sejc 
101.34 to 1689.0 ft/sec 

6 

±30° 

<P 

±90 rad 

q 

±20°/ sec 

r 

±10°/sec 

P 

±40°/sec 

O' 

-5° to +30° 

/3 

±20° 

h 

0 to 21 000 m 
0 to 63 000 ft 

®e 

+ 7° to -28° 


TABLE IV.- STANDARD DEVIATIONS OF CALCULATED STATES FROM 
MEASURED STATES AT CONVERGENCE 


“t» 

deg 

M 




State 




u 

w 

deg/sec 

deg 

a X’ 

g 

a Z’ 

g 

m/sec 

ft/sec 

m/sec 

ft/sec 

9.2 

0.71 

0.824 

2.7 

0.300 

0.98 

0.28 

0.17 

0.014 

0.028 

13 

.66 

1.98 

6.49 

.427 

1.4 

.21 

.11 

.007 

.033 

8.2 

.81 

.619 

2.03 

.372 

1.22 

.22 

.15 

.012 

.039 

12 

.78 

.573 

1.88 

.741 

2.43 

.62 

2.27 

.012 

.084 
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TABLE V.- PARAMETER VALUES EXTRACTED (STANDARD DEVIATIONS IN PARENTHESES) 


0.067 

0.0! 

2 x 10" 4 ) 

(±1.5 x 

0.40 

0.2! 

: 0.016) 

(±0.0C 

-0.53 

- 0.1 

0.0016) 

(±6.3 x 


t = 8° a t = 12° 
0.7 M = 0.7 M = 0.8 M = 0.8 


^Reference values are for an F-8C at 12.19 km (40 000 ft) and with the weight about 5 percent greater and 
the inertias about 15 percent greater than those of the subject aircraft. They have been transferred from a 
0.287c center of gravity to a 0.29c center of gravity and from stability axes to body axes. 




TABLE VI.- CORRELATION MATRICES 


t! 

C 

a 

a 

o 

Cn 


Run with a - 

= 9.2°, M = 0.71 


-0.058 

-0.028 

-0.0 

8 

1 

.18 

-.4 

8 

.18 

1 

-.0 

4 

-.47 

-.0022 

1 

0 

! 

• 

CO 

o 

.26 

.7 


Run with a ■ 

= 12°, M = 0.78 


-0.15 

-0.13 

0.1 


1 

.38 

-.7 


.38 

1 

-.2 


-.76 

-.27 

1 


-.08 

.17 

.5 


Run with a = 13°, M = 0.66 


-0.08 

-0.07 

1 

.29 

.29 

1 

-.57 

-.12 

-.25 

.16 

. 



Run with a = 8.2°, M= 0.81 
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TABLE VH.- PERIODS AND TIME TO DAMP TO HALF- AMPLITUDE 
FOR THE SHORT- PERIOD MODES OF THE RUNS MADE 


M 

a, 

deg 

t l/2,f’ 1 

sec 

t l/2,c» 

sec 

Pf, 

sec 

Pc, 

sec 

5c 

Sf 

0.71 

9.2 

1.1 

0.96 

2.20 

2.38 

0.260 

0.220 

.66 

13 

1.35 

1.20 

2.30 

2.31 

.208 

.192 

.81 

8.2 

1.20 

1.12 

2.70 

2.62 

.250 

.240 

.78 

12 

1.40 

1.36 

2.90 

2.75 

.218 

.217 

* .7 

9 


1.3 


2.4 

.200 



13 


1.4 


2.3 

.180 


*.8 

8 


1.2 


2.3 

.207 



12 


1.3 


2.4 

.200 



From reference 6. 
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Figure 2.- Comparison of measured data with time histories computed by using parameters 
given in table V for flight data taken at M = 0.71 and = 9.2°. 
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Figure 3.- Comparison of measured data with time histories computed by using parameters 
given in table V for flight data taken at M = 0.66 and a, = 13°. 
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Figure 4.- Comparison of measured data with time histories computed by using parameters 
given in table V for flight data taken at M - 0.81 and a t = 8.2°. 
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Figure 7.- Concluded 
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Figure 8.- Comparison of measured data not used during parameter extraction with 
time histories computed by using parameters of table V determined for flight 
conditions M = 0.66 and «i = 13°. 
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Figure 9.- Comparison of measured data not used during parameter extraction with 
time histories computed by using parameters of table V determined for flight 
conditions M = 0.81 and ou = 8.2°. 
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Figure 9.- Concluded. 
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